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Measurement of event shapes in deep 
inelastic scattering at HERA 

ZEUS Collaboration 


Abstract 

Inclusive event-shape variables have been measured in the current region of the Breit 
frame for neutral current deep inelastic ep scattering using an integrated luminosity of 
45.0 pb“^ collected with the ZEUS detector at HERA. The variables studied included 
thrust, jet broadening and invariant jet mass. The kinematic range covered was 10 < 
< 20480GeV2 and 6 x lO'^ < X < 0.6, where is the virtuality of the exchanged 
boson and x is the Bjorken variable. The Q dependence of the shape variables has been 
used in conjunction with NLO perturbative calculations and the Dokshitzer-Webber non- 
perturbative corrections (‘power corrections’) to investigate the validity of this approach. 
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1 Introduction 


The hadronic hnal states formed in e’''e“ annihilation and in deep inelastic scattering 
(DIS) can be characterised by a number of variables that describe the shape of the event. 
These variables are also defined at the parton level, where they are calculable using 
perturbative QCD (pQCD). A comparison of the hadron-level measurements with the 
parton-level calculations tests pQCD theory as well as QCD-based models. It is critical, 
however, that the non-perturbative effects due to hadronisation are correctly taken into 
account. 


A phenomenological determination of the hadronisation corrections can be obtained through 
the use of Monte Carlo (MC) models. However, an analytic method has been presented 
by Dokshitzer and Webber P] which allows the necessary corrections to be explicitly 
evaluated |2]. The mean value of a given shape variable is taken to be the sum of two 
parts, one of which is calculable perturbatively in QCD, while the other models the soft, 
non-perturbative contribution. The measured shape variable then depends on two experi¬ 
mentally determined constants, namely the QCD coupling parameter, as, and an effective 
non-perturbative coupling, c^. The non-perturbative contribution is an analytic expres¬ 
sion, known as a power correction P], which varies inversely with the hard interaction 
scale of the event. In annihilation, this scale is taken to be ^/s, the centre-of-mass 
energy of the incoming particles, while in DIS it is taken to be Q, the square root of the 
virtuality of the exchanged boson. The power corrections provide a potentially powerful 
tool in the experimental study of parton physics using perturbative QCD. They give a 
good description of event shapes in e’''e“ annihilation pilhllH] . A similar success was antici¬ 
pated in DIS, assuming the universality of quark fragmentation, and hence of event-shape 
properties. The hrst such studies have been reported by the HI Collaboration [7j. 


In annihilation, the event shapes may be evaluated in the laboratory frame. In order 
to study quark fragmentation in DIS, a frame that isolates the current-quark region of 
the event from the proton-remnant region is required, since only the current-quark region 
is of interest. A natural frame for this purpose is the Breit frame |H]. In this frame, the 
longitudinal axis is the direction of the incoming proton, and the current quark emerges 
in the opposite direction; the hnal-state particles are assigned to the current region if 
their longitudinal momentum component is negative, in which case they are interpreted 
as products of the hadronisation of the current quark. 


This paper presents measurements of event shapes in DIS. The validity of the power 
correction method is studied by examining whether the data can be correctly described 
by the theoretical expressions, and by checking the consistency of the values of and 
Oo obtained from hts to the different event-shape variables. Measurements are given of 
the mean values of the selected event-shape variables, evaluated in the kinematic range 
6 X 10“^ < X < 0.6, 10 < < 20 480 GeV^ and 0.04 < y < 0.95. Here x is the Bjorken 

variable and y = Q‘^/ sx. Dehnitions of the event shapes are given in the following section. 
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The Q dependence of the means of the event-shape variables is fitted to next-to-leading- 
order (NLO) estimates from pQCD pHITIl] . using the Dokshitzer-Webber power corrections, 
to determine values of ag and oq. 


2 Event-shape variables 


The event-shape variables studied here are thrust, jet broadening, the invariant jet mass 
and the C-parameter. 


Thrust measures the longitudinal collimation of a given hadronic system, while broadening 
measures the complementary aspect. These two parameters are specihed relative to a 
chosen axis, denoted by a unit vector n. Thus: 


T 


Ei \Pi • «l _ Ei bii 

Ei Ip* I Ei Ip* 


EJ 

Ip* 

; X 

n 

E*p± 

EJ 

|P*I 

EJ 

|p*l 


The sums in the formulae are taken over all particles in the chosen region of the event, 
namely the current region of the Breit frame. With n taken to be the virtual-photon 
direction, thrust and broadening are denoted by and respectively. Alternatively, 
both quantities may be measured with respect to the thrust axis, dehned as that along 
which the thrust is maximised by a suitable choice of n. In this case, the thrust and 
broadening are denoted by Tt and Bt- 


The normalised jet invariant mass, M, is dehned by 




(2E.B.)'' 


The C-parameter is given by 


C* — 3 (A1A2 -l- A2A3 -|- A3A1), 


( 1 ) 


the cyclic sum of the products of the eigenvalues of the linearized momentum tensor, 0"^, 
where 


0a/3 


Eii 

(pi“p 

■■1 

' / |p*l 

EJ 

Ip* 



Taking 6ij as the angle between two particles, Eq. 0 can be simplihed to 


CO 

lpil 

lpjl 

sin^ 9ij 

2(E, 

lpil)^ 
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As seen from the above equations, the shape parameters in the present study are nor¬ 
malised to the energy in the current hemisphere. With this normalisation, to ensure 
infra-red safety, it is necessary to exclude events in which the energy in the current hemi¬ 
sphere is less than a certain limit, Sum- Values of Eum of O.IQ and 0.25Q have been 
used. The primary analysis is based on event shapes calculated in the P-scheme, i.e. with 
particles taken to have zero mass after transformation to the Breit frame. The P-scheme, 
in which particle masses are assumed nn, was used as a cross-check. 

In the Born approximation, and Tt are unity. Consequently, the shape variables 
(1 — Ty) and (1 — Tt) are employed so that non-zero values at the parton level are a direct 
indicator of higher-order QCD effects. 


3 Detector description and event selection 


The data sample presented here was collected with the ZEUS detector during 1995-1997 
and corresponds to an integrated luminosity of 45.0±0.7 pb“^. During this period, HERA 
operated with protons of energy Ep = 820 GeV and positrons of energy Pe = 27.5 GeV. 

The ZEUS detector is described in detail elsewhere 112 ]. The main components used in 
the present analysis are the uranium-scintillator calorimeter (GAL) and the central 
tracking chamber (GTD) jHj, which are positioned in a 1.43 T solenoidal magnetic held. 
The GAL is divided into forward, barrel and rear sections,^ each of which is subdivided into 
cells whose energy deposits are read out independently. The relative energy resolution, 
as measured in test beams, is 18%/-\/P( GeV) and 35%/'yp/GeV7 for electrons and 
hadrons, respectively. The interaction vertex was measured using the GTD with a typical 
resolution of 0.4 (0.1) cm along (transverse to) the beam direction. 

The DIS kinematic variables and the four-vector of the virtual photon were reconstructed 
using the double-angle (DA) method ^3]. The following additional experimental quanti¬ 
ties are dehned: 

• 6 = SjPj(l — cos6'i), where Pj is the energy measured by a calorimeter cell, and 6i 
is the polar angle of the cell relative to the vertex position. The sum runs over all 
calorimeter cells; 

• 2/jb = EjPi(l — cos6i)/2Ee, where the sum runs over all calorimeter cells except those 
associated with the scattered positron; 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the ‘forward direction’, and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The laboratory polar 
angle, 0, is measured with respect to the proton beam direction. 
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• I/e = 1 — {E'^/2Ee){l — cos 6*e), where E'^ and 6^ are the energy and angle of the scattered 
positron. 

A pure sample of DIS events in a well-dehned kinematic region was selected as follows: 

• the events must pass the trigger, whose critical component in this case was a selection 
on a high-energy scattered positron identihed in the CAL; 

• a primary vertex must be reconstructed, with one or more tracks, satisfying —50 < 
Zvtx < 40 cm; 

• a well-identihed scattered positron, which was found by means of a neural-network 
procedure uni, must satisfy E'^ > 10 GeV; 

• the measured impact position of the scattered positron on the face of the rear calorime¬ 
ter must be outside a square of 16 x 16cm^ centred on the beampipe; 

• 35 < 5 < 60 GeV, to remove photoproduction events where the scattered positron is 
lost down the beampipe, and also to reduce the effects of initial-state radiation; 

• 2/jb > 0.04, to ensure a well-measured hadronic system; 

• I/e < 0.95, to reduce background contamination from photoproduction HH; 

• > 10 GeV^ and x > 6 x lO”'^. 

Both track and calorimeter information were used to determine the event shapes. Galorime- 
ter cells were hrst grouped to form clusters and these clusters were then associated 
with tracks, where possible, to form energy-flow objects, EFO’s, associated with the 
hadrons formed in the interaction HSj. Both tracks and clusters were required to have 
20° < 9 < 160° and transverse momentum > 150 MeV. At this stage, EFO’s with 
tracks were assigned the mass of the pion, while those without tracks, corresponding 
mainly to photons from vr® decays, were assigned zero mass. Each accepted EFO was 
then transformed to the Breit frame, where it was assigned to the current region if its 
longitudinal momentum was negative. Subsequently, the masses were assigned according 
to the P- or E-scheme. A total of 321 000 events resulted from these selections and were 
used in the analysis. 


4 QCD models and event simulation 

Monte Garlo event simulation was used to correct the data for acceptance and resolution 
effects. The detector simulation was performed with the GEANT 3.13 program na. 
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Neutral current DIS events were generated using the DJANGOH 1.1 package [201; com¬ 
bining the LEPTO 6.5.1 |2I] generator with the HERACLES 4.6.1 program [221; which 
incorporates hrst-order electroweak corrections. The parton cascade was modelled with 
the colour-dipole model (CDM), using the ARIADNE 4.08 program. In this model, 

coherence effects are implicitly included in the formalism of the parton cascade. The Lund 
string-fragmentation model |2S] is used for the hadronisation phase, as implemented in 
JETSET 7.4 [20]. 

Further samples were generated with the HERWIG 5.9 program izzi, which does not 
apply electroweak radiative corrections. The coherence effects in the hnal-state cascade 
are included by angular ordering of successive parton emissions, and a clustering model 
is used for the hadronisation |2H]- Events were also generated using the MEPS option of 
LEPTO within DJANGOH, which subsequently uses a parton-showering model similar 
to HERWIG. To achieve agreement with the data, a diffractive component of 14% of the 
DIS events [221 was required and was simulated using RAPGAP 2.08 EDI- 

For ARIADNE, the default parameters were used. The LEPTO simulation was run with 
soft-colour interactions turned off, and HERWIG was retuned^ to give closer agreement 
with the measured shape variables at low Q; the GTEQ4D EH parameterisations of the 
proton parton distribution functions (PDF) were taken. The Monte Garlo event samples 
were passed through reconstruction and selection procedures identical to those for the 
data. 


5 Data correction 

The event shapes were evaluated for event samples in selected bins of x and Q^. The choice 
of the bin sizes EH was motivated by the need to have good statistics while keeping the 
migrations, both between bins, and from the current to the target region within each bin, 
to a minimum. The kinematic bin boundaries are listed in Table H 

In each (x, Q‘^) bin, the ARIADNE MG was used to investigate the event acceptance and 
the acceptance in each bin of the event-shape variable. The acceptance was dehned at 
the hadron level as the ratio of the number of reconstructed and selected events to the 
number of generated events in the given bin. The event acceptance exceeded 70% for all 
bins, while the event-shape acceptance was less than 70% only at the extremes of the 
range and at low y. 

Agreement was found between the uncorrected data and the predictions of ARIADNE 
throughout the entire kinematic range of each event-shape variable, thus conhrming its 
suitability for the purposes of correcting the data. Good agreement with ARIADNE 
was also found for the energy-flow [22] and charged-track distributions |22| studied in 

^ The parameter PSPLT was set equal to 1.8; otherwise default parameters were used. 
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previous analyses. The data were also compared with the HERWIG predictions; here the 
agreement with data was satisfactory but less good than for ARIADNE. 

The correction factors for the means of the shape variables were evaluated as the ratios 
of the generated to the observed values of the mean in each [x, Q^) bin. These correction 
factors, which were used for the subsequent analysis, lie within the range 0.75 - 1.12 and 
are typically within ±5% of unity. As a check, the calculation was repeated correcting 
the individual bins using the acceptances as described above. 

The generated distributions include the products of strong and electromagnetic decays, 
together with and A decays, but exclude the decay products of weakly decaying 
particles with lifetime greater than 3 x 10“^°s. The correction procedure accounts for 
event migration between [x, intervals, QED radiative effects, EFO-reconstruction 
efficiency and energy resolution, acceptances in pt and 6, EFO migration between the 
current and target regions, and the decay products of Kg and A decays that were assigned 
to the primary vertex. 


6 Systematic checks 

The systematic uncertainties in the measurement can be divided into three types, due to 
the MC model used, to the event reconstruction and selection, and to the EFO recon¬ 
struction. The systematic checks were as follows: 

• the data were corrected using a different hadronisation and parton-shower model, 
namely HERWIG or LEPTO, in place of ARIADNE; 

• the cut on yj^ was increased to 0.05; 

• the cut on 6 was changed to 40 < 5 < 60 GeV; this harder cut estimates any residual 
uncertainties in the photoproduction background; 

• the double-angle kinematics were recalculated after removing GAL deposits due to 
backwardly scattered particles (albedo) from the material close to the proton beampipe 

ESES]; 

• the measured energies of clusters in the calorimeter were varied by ±3%, ±1% and 
±2% for the forward, barrel and rear GAL sections, respectively, corresponding to the 
uncertainties on the associated energy scales; 

• the EFO cuts at 6^ = [20°, 160°] and pt > 150 MeV were tightened to 0 = [25°, 155°] 
and Pt > 200 MeV; the cuts were also removed; 

The largest systematic uncertainty was due to using HERWIG as the hadronisation model. 
The LEPTO model produced smaller changes. The other systematics were typically at 
the level of 1%, smaller than or similar to the statistical errors. 
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7 Power corrections 


Next-to-leading-order QCD calculations of the shape parameters have been made using 
the programs DISASTER++ and DISENT ^U], which give parton-level distributions. 
Both programs used the CTEQ4A PDFs. To determine the theoretical ag dependence 
of the variables, high-statistics event samples were generated for each of hve values of 
as{Mz) [2^. For both NLO calculations, the mean value of each shape variable was 
found to be linearly dependent on a* in the appropriate range. For each bin, therefore, 
the calculated value of the shape variable may be used to estimate the value of ag at the 
chosen reference scale, namely the mass of the boson. 

The relationship used to calculate ag{Q) in terms of ag{Mz) is [211 


ag{Q) = 


ag{M^ 


l + a.(Af.)L(”) -2- 


where, for the 2-loop form. 


L(2) 


Q 


+ A) log ^ ^ 


27r Svr^ 




with 


^ ^ 306 - 38Nf 

f3o = ll- -Nf and /3i =- - -. 


Here, Nf is the number of active quark flavours at the scale Mz, taken to be hve. 


Before comparison with experimental data, the calculated values of the shape parameters 
require correction for the effects of hadronisation. Dokshitzer and Webber calculated 
power corrections to the event-shape variables in e’''e“ annihilation, assuming an infrared- 
regular behaviour of the effective coupling, aes mq. The technique was subsequently 
applied to the case of DIS [3] and has been used here. 


In this approach, a constant, oq, is introduced, which is taken to be independent of the 
shape variable. This constant is dehned as the hrst moment of the effective strong coupling 
below the scale p/ and is given by: 


_ 1 /"w 

h/ Jo 


where the variable /x/ is the lower limit for the perturbative approach to be valid. This is 
taken to be 2 GeV, in common with previous analyses jHIHlinilZ! • 

The theoretical prediction for a mean event-shape variable, denoted by {V), is then given 
by 

(V) = {v}^^o + {r>p„w. (-3) 
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where (^)nlo calculated perturbatively, and is the power correction, 

generalized power correction is given by 


{V) 


pow 


dv 


AMAi 

nQ 


The 


(- 3 ) 


This has a 1/Q dependence with a calculable coefficient, ay- The variable M. is the ‘Milan 
factor’ of value 1.49 isn, which takes into account two-loop corrections; it has a relative 
uncertainty of about ±20%, due to three- and higher-loop effects. The term Ai is given 
by: 


= —/i/ 
TT 


tto - asinn) - ^flogf — 


27r 


V i^i 


K 


Po 


- 3 ) 


where Cp = K = 


^ \— |-/V/ and Nf is taken to be five. The central analysis was 

performed with the renormalisation scale /i^ set to Q. The dependence of the perturbative 
predictions on this scale was studied by introducing a variable parameter, such that 
= x^Q. The NLO calculations, but not the power corrections, depend also on the 
factorisation scale, which is varied from its central value Q by introducing a similar 
parameter, xp- 


In Eq. ((Zl), the low-energy contribution to the mean shape variable is determined by 
Oo, while the remaining terms subtract out the integral, up to p/, of the perturbative 
expression for the shape average. Above this limit, the perturbative expression is taken 
to be applicable. 

The values of ay for (1 — T^), (1 — Ty), C and are respectively 2, 2, 1 and Svr. For 
i?y, the form |2HI 

ay = — _ ± 0.75 — ^ ± (-3) 

2-^ 2C'i?acMw(Q) ^ 

has been used, where ? 7 o = —0.614, Q = and the physical coupling acMW is related 

to the standard ajjg by 

ttcMW = • 

In Eq. ©, Pq was calculated with Nf = 3, which is appropriate for fragmentation. In 
the corresponding expression for Bp, the factors 2Cp and 12Cp were replaced hy Cp and 
QCp, respectively. 


It has recently been found IHH] that for 5^, there is an additional x-dependent term in 
Eq. ((7j). The importance of this term in the present analysis has been found to be small; 
for further comments, see below. 


8 Results 

In this section, the observed Q dependences of the mean values of the event-shape pa¬ 
rameters are presented and compared with the expectations from theory. To study the 









sensitivity to x, results for Q'^ < 320 GeV^ are calculated separately in the high and low 
X ranges as well as for the full x range. The values of Sum used were 0.25Q, as recom¬ 
mended goi to ensure convergence of the perturbative series involving \n{Q /2Siim) , and 
O.IQ HI. The mean values of the event-shape parameters for Sum = 0.25Q are listed in 
Table El together with their statistical and systematic uncertainties. 

Figure^ shows the corrected mean values of the event shapes as a function of Q, using 
^lim = O.IQ, together with the HI measurements [7j, with which there is good agreement. 
The mean values fall with Q at the higher Q values, and at lower Q show an x dependence 
at hxed Q which is most pronounced for the variables measured with respect to the photon 
axis. Both here and with Sum = 0.25Q (not shown), good agreement is found with the 
predictions of ARIADNE over the entire Q range. The agreement with HERWIG is less 
good, particularly for the variables (1—Tr), Bt, and C at low Q. Although not used in 
the subsequent hts, the data for < 80 GeV^ are included for purposes of completeness, 
and to allow comparison with more extended theoretical hts if they become available. 

The data were htted to the sum of an NLO term, obtained from DISASTER-I--I- or 
DISENT, plus a power correction as described above. The theoretical calculation neglects 
terms of the order of 1/Q^ p]. Gonsequently, the analyses have been conhned to the region 
Q > 9 GeV, i.e. bins 7 - 16 in Tabled With fij hxed at 2 GeV, there are two parameters 
in Eq. 0, as{Mz) and ao, that can be varied to obtain the best agreement between 
calculation and data. 

Two types of ht were studied: the ohset method and the Hessian method. The ohset 
method |lllll2j uses a dehned using a diagonal error matrix with errors given by the 
statistical errors on the data points combined in quadrature with the uncorrelated errors 
arising from the limited statistics of the DISENT and DISASTER-I--I- calculations. The 
Hessian method HSUH uses a full error matrix which includes correlated oh-diagonal 
terms due to the systematic uncertainties. The hts obtained using the ohset method 
are shown in Fig. El For the Hessian method, the four major systematics, namely those 
associated with S, the tightening and relaxation of the EFO angular and pt cuts, and the 
use of HERWIG, were included as oh-diagonal terms in the error matrix. As expected, 
the is reduced with the Hessian method, and the htted error, which includes the 
systematic contribution, is approximately a factor of two smaller than the systematic 
uncertainties estimated from the ohset hts. For the variables (1 — Tt), Bt-, M"^ and G, 
the values of from the two ht methods agree within the statistical uncertainties. For 
the tts values from (1 — T,y) and B,y, and for all evaluations of Oq, the two ht methods 
give results that in general diher by more than the statistical errors; they do however 
agree within the ohset systematic uncertainty. The diherence in the results of the hts 
originates primarily from the use of HERWIG in place of ARIADNE. The Hessian method 
relies upon Bayesian priors, specihcally the assumption of Gaussian distributions, for the 
systematic uncertainties. There is no reason to believe that this assumption is correct for 
the dominating fragmentation systematic. Gonsequently, the results presented here are 
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based on the offset method with its conservative estimate of systematic effects. 

The DISASTER++ calculation gives predictions in closer agreement with analytic calcu¬ 
lations |4()[l44j than does DISENT, which is believed to contain errors |33]. The hts using 
DISASTER-I--I- for the NLO term are therefore taken to give the more reliable estimates 
of as and oq. The DISENT-based analysis provides a check and facilitates comparison 
with HI, whose analysis used this program. Reasonable hts are obtained; those based on 
DISASTER-!--!-, using Eum = 0.25Q, are shown in Fig. |2l while those for DISENT (not 
shown) are similar. The htted power-correction term is substantial except for the variables 
that are based on the virtual-photon axis. The results of the hts using DISASTER-!--!- 
and DISENT are compared in Fig. El where it can be seen that, with the exception of 
as is the same within statistical errors for the two NLO calculations. For all variables, Oq 
determined using DISASTER-!--!- is smaller than when using DISENT. 

As seen in Fig. El the data have a signihcant x-dependence at a given Q; consequently, 
hts have been made using the high-x and low-x selections, as well as to the full set of used 
bins. Results from the DISASTER-!--!- hts are shown in the contour plots of Fig. EJ 

To allow a direct comparison with the HI results [Z], Fig. El shows the determinations of a* 
and Oo using DISENT and Sum = O.IQ; as expected from the agreement of the measured 
data points (Fig. (H), the agreement in the htted as and oq values is, in general, good. In 
particular, a low value of as for is conhrmed. The inhuence of the Sum selection is 
also illustrated in Fig. El which shows, for DISENT, that the diherent selections lead to 
Oq values that agree within about two standard deviations. In contrast, appears more 
sensitive to Sum for several of the variables. In general, it is found that Sum = 0.25Q 
gives a smaller variation of the htted a* and oq with x than is found for Sum = O.IQ (not 
shown). Given this, together with the theoretical preference for the higher Sum (101 , the 
central analysis is based on the data evaluated with Sum = 0.25Q. 

The experimental systematic uncertainties on as and aH were estimated by repeating 
the hts with the systematic variations described in Section El The largest ehect resulted 
from correcting the data with HERWIG instead of ARIADNE. The use of HERWIG gave 
a systematic increase in as and a systematic decrease in oq; these shifts are possibly 
attributable, respectively, to the use of parton showers rather than the colour-dipole 
model, and to the diherent fragmentation schemes used in the models. Also, HERWIG 
does not contain electroweak terms. The variables (1—T..y) and B.y are in addition sensitive 
to the method of reconstructing the kinematic variables, owing to their dependence on 
the photon direction. 

Figure El summarises the ag and values from hts using the statistical errors, in order 
to indicate the degree of agreement between the diherent measurements. The systematic 
uncertainties on the data introduce highly correlated ehects on the results obtained from 
the diherent shape variables, and so are not included here. The inconsistency which is 
evident between the diherent determinations is discussed below. 
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To estimate the theoretical uncertainties, the fragmentation and renormalisation scales 
were varied by a factor of two, and studies were made of the effects of changes to /i/ and to 
the Milan factor. To give an indication of the uncertainties due to mass effects, the data 
were reanalysed using the i?-scheme m It was found that depends strongly on /i/, 
decreasing as fij is increased. If the model were robust, as should have little dependence 
on /x/. 

The £t results including experimental and theoretical systematic uncertainties are col¬ 
lected in Tables 121 and The dominant uncertainty comes from the variation of the 
renormalisation scale. The renormalisation-scale uncertainty quoted here follows the pro¬ 
cedure employed in the e+e“ studies iHEiini using Eq. ©• If, following HI [7j, Q in 
Eq. (jTj) were replaced by p/j, the theoretical uncertainty due to the variation would 
be approximately a factor of two larger. The influence of the a:-dependent term in Eq. 
© was examined using the CTEQ5M proton structure jlH]- While an improved fit to 
the data was obtained, the changes to the resulting as and oq values were less than their 
statistical uncertainties. 


9 Discussion 


From the results using DISASTER-|--|-, the following features are observed: 

• (1 — Ty) requires a smaller hadronisation correction than (1 — Tt) (Fig. |21), contrary to 
the theoretical expectation that the correction should be similar. This is responsible 
for the signihcantly different oq values for the two thrust variables; 

• (1 — Ty) shows a larger a:-dependence than (1 — Tt)'-, the and C variables, whose 
dehnition does not depend on a choice of axis, show a small x dependence (Fig. 12]); 

• a residual x dependence in the htted oq value obtained from (1 — Tt), (1 — Ty) and 

but not from Bt, B,y and C can be seen in Fig. 01 However, the Oq values are 
also consistent with a similar small x dependence in all four variables; 

• the as values from all the variables except (1 — Ty) show no significant x dependence. 

There is an inconsistency between the Oq values determined from (1 — Ty) and from the 
other variables. As noted earlier, there is an x dependence in the data from (1 — Ty) that 
is not well described by the present model (Fig. 0|). Consequently, the htted values of as 
and Oo for this variable are unlikely to be meaningful. The anomalous value of a* may 
also be due to the fact that DISASTER-I--I- does not take full account of initial-state gluon 
radiation or other effects related to the target remnant, which may affect the direction of 
the current-region system. 
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For the variables (1 — T^), and C, the htted values are consistent within the 

statistical uncertainties. The variable Bx gives an value that differs from the other 
determinations, although its oq value agrees within ±10% with those from the other four 
variables. The inconsistency cannot be due to experimental or theoretical systematic 
uncertainties, including scale uncertainties, since these act in the same direction for all 
the variables. It may be taken to indicate that Bj- has a greater sensitivity to higher-order 
corrections than the other variables. 

A comparison with other measurements is of interest. With the exception of (1 — Ty), 
the present results for as and oq are consistent with those from data within the 

substantial theoretical uncertainties. Using data from a variety of experiments, 

Movilla Fernandez et ah found good agreement between the means of the event-shape 
variables as a function of a/s; in contrast with the observations of this paper, they obtained 
as values that were consistent within statistical errors for all the variables studied nnai. 
The Cio values were likewise mutually compatible, with the possible exception of that from 
M^. An overall consistency was claimed which confirmed the validity of the model in the 
context and thus enabled an overall experimental value for ag to be given. 

In summary, the power correction method applied in DIS gives consistent values for as 
for the event-shape variables (1 — Tp), B^, and C. The Oq values for these variables 
agree to within ±10%, which is consistent with the precision claimed for the model. The 
variables (1—Ty) and Bx give, respectively, Oq and a* values that are inconsistent with the 
other determinations. It must be concluded, therefore, that the power-correction model 
does not consistently describe all the shape variables in DIS. Consequently, no average as 
or oo values are quoted. 


10 Summary 

A measurement has been made of the mean values of the event-shape variables thrust (T), 
broadening (B), normalised jet mass (M^) and C-parameter, using the ZEUS detector 
at HERA. The variables T and B were determined relative to the virtual photon axis 
and the thrust axis. The events were analysed in the Breit frame for the kinematic range 
6 X10“"^ < X < 0.6, 10 < < 20 480 GeV^ and 0.04 < y < 0.95. The data are successfully 

described by the ARIADNE Monte Carlo model. 

The Q dependences of the mean event shapes have been fitted to NLO calculations 
from perturbative QCD with the DISASTER±± and DISENT programs together with 
the Dokshitzer-Webber non-perturbative power corrections, with the aim of determining 
as{Mz) and cF)(h/)- Such a model should give values of ag and Oq that are independent 
of the shape variable. Neither DISASTER±± nor DISENT fulfils these requirements for 
all variables. 

Using DISASTER±±, consistent values of ag are obtained for the shape variables (1—T^), 
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5^, and C, with values that agree to within ±10%. The Oq value from (1 — T^) 
and the ag value from are in disagreement with the other determinations. With the 
exception of (1 — Ty), the present values are consistent with those measured in e’''e“ 
annihilation, to within the theoretical uncertainties. There is consistency with the results 
from HI. 

The power correction method provides a successful description of the data for all event- 
shape variables studied. Nevertheless, the lack of consistency of the ag and oq determina¬ 
tions obtained in deep inelastic scattering, together with the dependence of the results on 
Bjorken x, suggest the importance of higher-order processes that are not yet included in 
the model employed in this analysis. These effects must be understood before a reliable 
value of ag can be quoted using the present method. 
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Bin 

g' (GeV^) 

X 

Bin 

g' (GeV^) 

X 

1 

10 - 20 

0.0006 - 0.0012 

9 

160 - 320 

0.0024 - 0.010 

2 

10 - 20 

0.0012 - 0.0024 

10 

160 - 320 

0.01 - 0.05 

3 

20 - 40 

0.0012 - 0.0024 

11 

320 - 640 

0.01 - 0.05 

4 

20 - 40 

0.0024 - 0.010 

12 

640 - 1280 

0.01 - 0.05 

5 

40 - 80 

0.0012 - 0.0024 

13 

1280 - 2560 

0.025 - 0.150 

6 

40 - 80 

0.0024 - 0.010 

14 

2560 - 5120 

0.05 - 0.25 

7 

80 - 160 

0.0024 - 0.010 

15 

5120 - 10240 

0.06 - 0.40 

8 

80 - 160 

0.01 - 0.050 

16 

10240 - 20480 

0.10 - 0.60 


Table 1: The kinematic boundaries of the bins in x and Q^. The power-correction 
fits use bins 7 - 16, apart from the fits denoted ‘high x’ and ‘low x\ which omit 
bins 7, 9 and bins 8, 10, respectively. 
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Bin 

1 — JT 

Bx 

m2 

C 

1 -Ty 


1 

0.1595 ± 0.0010 ± 0.0046 

0.2172 ± 0.0009 ± 0.0060 

0.0967 ± 0.0006 ± 0.0039 

0.540 ± 0.003 ± 0.014 

0.508 ± 0.002 ± 0.009 

0.3932 ± 0.0010 ± 0.0066 

2 

0.1609 ± 0.0010 ± 0.0041 

0.2185 ± 0.0009 ± 0.0059 

0.0965 ± 0.0006 ± 0.0042 

0.544 ± 0.003 ± 0.013 

0.499 ± 0.002 ± 0.018 

0.3892 ± 0.0010 ± 0.0108 

3 

0.1672 ± 0.0011 ± 0.0078 

0.2223 ± 0.0010 ± 0.0094 

0.0924 ± 0.0007 ± 0.0054 

0.565 ± 0.003 ± 0.023 

0.458 ± 0.003 ± 0.003 

0.3731 ± 0.0011 ± 0.0013 

4 

0.1708 ± 0.0009 ± 0.0102 

0.2248 ± 0.0008 ± 0.0123 

0.0941 ± 0.0005 ± 0.0079 

0.575 ± 0.002 ± 0.031 

0.437 ± 0.002 ± 0.004 

0.3644 ± 0.0009 ± 0.0044 

5 

0.1620 ± 0.0020 ± 0.0070 

0.2155 ± 0.0014 ± 0.0056 

0.0871 ± 0.0009 ± 0.0032 

0.556 ± 0.004 ± 0.015 

0.426 ± 0.004 ± 0.003 

0.3580 ± 0.0020 ± 0.0020 

6 

0.1679 ± 0.0010 ± 0.0114 

0.2204 ± 0.0009 ± 0.0122 

0.0880 ± 0.0006 ± 0.0060 

0.572 ± 0.003 ± 0.033 

0.3968 ± 0.0024 ± 0.0015 

0.3452 ± 0.0011 ± 0.0012 

7 

0.1500 ± 0.0004 ± 0.0050 

0.2040 ± 0.0004 ± 0.0049 

0.0791 ± 0.0002 ± 0.0030 

0.5314 ± 0.0010 ± 0.0143 

0.3547 ± 0.0010 ± 0.0031 

0.3238 ± 0.0005 ± 0.0016 

8 

0.1536 ± 0.0005 ± 0.0110 

0.2061 ± 0.0004 ± 0.0097 

0.0808 ± 0.0003 ± 0.0084 

0.5403 ± 0.0013 ± 0.0300 

0.3149 ± 0.0012 ± 0.0190 

0.3031 ± 0.0006 ± 0.0096 

9 

0.1322 ± 0.0007 ± 0.0013 

0.1866 ± 0.0006 ± 0.0020 

0.0703 ± 0.0004 ± 0.0008 

0.483 ± 0.002 ± 0.006 

0.3130 ± 0.0020 ± 0.0050 

0.3009 ± 0.0009 ± 0.0029 

10 

0.1347 ± 0.0006 ± 0.0054 

0.1875 ± 0.0005 ± 0.0053 

0.0718 ± 0.0003 ± 0.0051 

0.489 ± 0.002 ± 0.017 

0.2742 ± 0.0013 ± 0.0071 

0.2789 ± 0.0007 ± 0.0045 

11 

0.1150 ± 0.0008 ± 0.0027 

0.1676 ± 0.0007 ± 0.0033 

0.0628 ± 0.0004 ± 0.0032 

0.431 ± 0.002 ± 0.010 

0.242 ± 0.002 ± 0.005 

0.2574 ± 0.0011 ± 0.0026 

12 

0.0982 ± 0.0012 ± 0.0017 

0.1489 ± 0.0011 ± 0.0024 

0.0543 ± 0.0007 ± 0.0017 

0.375 ± 0.003 ± 0.006 

0.210 ± 0.003 ± 0.004 

0.235 ± 0.002 ± 0.002 

13 

0.086 ± 0.002 ± 0.002 

0.133 ± 0.002 ± 0.002 

0.0487 ± 0.0010 ± 0.0021 

0.330 ± 0.005 ± 0.006 

0.168 ± 0.004 ± 0.008 

0.203 ± 0.003 ± 0.003 

14 

0.078 ± 0.003 ± 0.004 

0.122 ± 0.003 ± 0.005 

0.043 ± 0.002 ± 0.002 

0.295 ± 0.009 ± 0.011 

0.120 ± 0.007 ± 0.013 

0.167 ± 0.005 ± 0.006 

15 

0.069 ± 0.005 ± 0.004 

0.112 ± 0.006 ± 0.004 

0.039 ± 0.003 ± 0.003 

0.260 ± 0.015 ± 0.012 

0.125 ± 0.012 ± 0.015 

0.160 ± 0.009 ± 0.012 

16 

0.062 ± 0.012 ± 0.015 

0.104 ± 0.013 ± 0.017 

0.030 ± 0.007 ± 0.009 

0.233 ± 0.035 ± 0.035 

0.058 ± 0.022 ± 0.031 

0.110 ± 0.018 ± 0.018 


Table 2: Mean event-shape variables in the bins dehned in Table ^ The hrst 
uncertainty is statistical and the second is systematic. 












Variable 

1 — Tp 

Bp 


c 

l-T^ 



0.1258 

0.1159 

0.1271 

0.1274 

0.1354 

0.1270 

stat. error 

±0.0013 

±0.0013 

±0.0016 

±0.0010 

±0.0028 

±0.0026 

stat.+sys. unc. 

±0.0040 

±0.0026 

±0.0040 

±0.0021 

±0.0132 

±0.0102 

X^/dof 

2.8 

1.3 

2.1 

1.5 

2.5 

1.8 

correlation 

-0.25 

-0.80 

-0.60 

0.26 

-0.93 

-0.84 

xp = 0.5 

-0.0007 

-0.0001 

-0.0004 

-0.0007 

±0.0088 

±0.0025 

xp = 2.0 

±0.0009 

±0.0008 

±0.0003 

±0.0007 

-0.0008 

±0.0036 

x^ = 0.5 

-0.0068 

-0.0067 

-0.0080 

-0.0067 

-0.0088 

-0.0222 

x^ = 2.0 

±0.0083 

±0.0081 

±0.0090 

±0.0082 

±0.0084 

±0.0053 

M = 1.19 

±0.0025 

±0.0018 

±0.0024 

±0.0029 

±0.0032 

±0.0014 

M = 1.79 

-0.0021 

-0.0017 

-0.0022 

-0.0025 

-0.0027 

-0.0012 

H-i = 1 GeV 

±0.0054 

±0.0042 

±0.0053 

±0.0063 

±0.0069 

±0.0029 

/i/ = 4 GeV 

-0.0059 

-0.0047 

-0.0061 

-0.0068 

-0.0075 

-0.0039 

E-scheme 

±0.0040 

±0.0028 

±0.0029 

±0.0030 

±0.0025 

±0.0010 

Total 

+0.0117 

+ 0.0101 

+ 0.0118 

+ 0.0114 

+ 0.0197 

+ 0.0128 

uncertainty 

-0.0101 

-0.0088 

-0.0111 

-0.0101 

-0.0178 

-0.0248 


Table 3: Fitted results for as{Mz) using the NLO prediction from DISASTER++ 
and Slim = 0.25Q. The quoted is that from the offset-method £t using sta¬ 
tistical uncertainties and DISASTER-I--I-. The third line is the total experimental 
uncertainty from the statistical and experimental systematic uncertainties added in 
quadrature. The quoted is from the offset-method fit using statistical uncertain¬ 
ties and DISASTER-I--I-. The fifth row gives the correlation coefficients between 
the fitted values of as{Mz) and oq (see next Table). The xp, x^, M, and fij rows 
give the theoretical systematic uncertainties due to variations on the fragmentation 
and renormalisation scales, the Milan factor and the lower limit for the perturba¬ 
tive calculation, respectively; the xp and x^ values denote factors by which the 
respective scale values are varied. The systematic effect of using the E-scheme 
rather than the P-scheme is given in the final row of systematic uncertainties. The 
total uncertainty is the total experimental uncertainty added in quadrature with 
the theoretical uncertainties. 
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Variable 

1 — Tp 

Bp 

M2 

C 

1 

B.y 

oh {2 GeV) 

0.4843 

0.4566 

0.4440 

0.4274 

0.3286 

0.4593 

stat. error 

±0.0020 

±0.0041 

±0.0030 

±0.0017 

±0.0187 

±0.0171 

stat.+sys. error 

±0.0264 

±0.0139 

±0.0439 

±0.0144 

±0.0993 

±0.0815 

xp = 0.5 

±0.0030 

±0.0019 

±0.0006 

±0.0018 

-0.1173 

-0.0491 

xp = 2.0 

-0.0017 

-0.0034 

±0.0002 

-0.0006 

±0.0421 

-0.0128 

x^ = 0.5 

±0.007 

±0.056 

±0.016 

±0.007 

±0.009 

±0.256 

x^ = 2.0 

-0.002 

-0.033 

-0.007 

-0.003 

±0.029 

-0.135 

M = 1.19 

±0.0363 

±0.0390 

±0.0223 

±0.0235 

-0.0402 

±0.0250 

M = 1.79 

-0.0282 

-0.0296 

-0.0192 

-0.0200 

±0.0187 

-0.0197 

E-scheme 

±0.0163 

±0.0101 

±0.0134 

±0.0127 

±0.0052 

±0.0103 

Total 

+0.0483 

+ 0.0706 

+0.0535 

+ 0.0312 

+ 0.1137 

+ 0.2700 

uncertainty 

-0.0387 

-0.0467 

-0.0484 

-0.0248 

-0.1589 

-0.1664 


Table 4: Fitted results for defined at fcj = 2 GeV, using the NLO prediction 
from DISASTER++ and Sum = 0.25Q. Other definitions are given in the caption 
to Table El 
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Figure 1: Comparison of the mean event-shape parameters (solid points) with 
ARIADNE and HERWIG predictions. The plots refer to a) thrust with respect to 
the thrust axis, b) jet broadening with respect to the thrust axis, c) invariant jet- 
mass squared, d) C-parameter, e) thrust with respect to the virtual photon axis, 
f) broadening with respect to the virtual photon axis. The data were corrected 
using ARIADNE. The inner error bars are statistical; the outer are statistical plus 
systematic added in quadrature. The open squares are HI data, with statistical 
uncertainties which are in most cases covered by the symbol. In e) and f), where 
the X variation is biggest, smoothed curves are drawn through the MC points for 
all X > 320 GeV^) and low x < 320 GeV^), and separately through the 
high-a; MG points < 320 GeV^). A value of Eum = O.IQ was used. 
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Figure 2: Fits to mean values of the shape variables versus Q, with Sum = 
0.25Q. Plots a) - f) are dehned in the caption to Fig. ^ The lines join the £t 
values at the Q values of the data: the solid line is the htted NLO prediction 
from DISASTER++ plus the power correction, while the dashed line is the htted 
DISASTER++ contribution alone. ‘High x’ and ‘low x’ refer to the subdivisions 
as dehned in Table 1; ‘all x’ refers to the points with Q > 20 GeV, which are not 
subdivided in x. 
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Figure 3: Contour plots for the parameters as{Mz) and oq fitted to the mean 
values of thrust and broadening measured with respect to the photon axis, jet-mass 
squared, C-parameter, and thrust and broadening measured with respect to the 
thrust axis. Results are shown for fits, using all data points for > SOGeV^, 
based on DISASTER-I--I- and DISENT with Sum = 0.25Q. The contours show the 
one-standard-deviation limits determined using statistical uncertainties only. 
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Figure 4: Contour plots for as{Mz) and Oq fitted to the mean values of the 
event-shape variables. The fits are based on DISASTER-I--I-, with Sum = 0.25Q. 
The contours show the one-standard-deviation limits determined using statistical 
uncertainties only. The high-x and low-x selections are as defined in Table 1, while 
‘all data’ uses all measured bins for > 80 GeV^. 
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Figure 5: Contour plots for as{Mz) and oq fitted to the mean values of the event- 
shape variables. The hts are based on DISENT, using all data for > SOGeV^, 
with energy cuts Sum = O.IQ (open squares) and 0.25Q (hlled circles). The contours 
show the one-standard-deviation limits determined using statistical uncertainties 
only. The rectangles enclose the associated pairs of points to guide the eye. 
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Figure 6: Contour plots for as{Mz) and oq fitted to the mean values of the 
event-shape variables. The fits are based on DISASTER-I--I-, with Sum = 0.25Q. 
The contours show the one-standard-deviation limits determined using statistical 
uncertainties only. For further comments, see text; the full systematic errors, which 
are strongly correlated between the different variables, are given in Tables 3 and 4. 
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